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Abstract 

The magnetocrystalline anisotropy energy E an i s for a monolayer of Fe and Ni 
is determined using a fully convergent tight-binding calculation including s-d 
hybridization. The spin-orbit interaction A so is treated non-perturbatively. 
Remarkably, we find E an is oc A^ and important contributions to E an i s due 
to the lifting of degeneracies near the Fermi-level. This is supported by the 
calculated decrease of the anisotropy energy with increasing temperature on 
a scale of several hundred K. Our results clarify the present debate on the 
origin of E anis . 
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Despite increasing theoretical effort, the origin of magnetic anisotropy in ferromagnetic 
transition metals has not been clearly identified. Actually, present theoretical analyses are 
controversial It remains to determine clearly how the magnetocrystalline anisotropy 

energy E anis depends on the spin-orbit coupling (SOC) strength A so and whether the level 
crossings at the Fermi-level Ep play an important role for E an i S or alternatively whether 
they have to be excluded as stated by Wang et al. 0. It is important to present calculations 
without using the state-tracking method by Wang et al. |1[|, the validity of which was disputed 
by Daalderop et al. 0. Due to increased recent experimental activity on magnetic anisotropy 
at surfaces, interfaces and thin films H|J, there is need of a theory without too many 
restrictive assumptions in order to get a systematic understanding of E anis . 

In this paper, we present a theory treating A so non-perturbatively and determining the 
electronic bandstructure within the combined interpolation scheme ||. We investigate a 
simple quadratic Fe and Ni monolayer epitaxially grown on the Cu(OOl) surface and neglect 
further interactions with the substrate. Our calculation permits us to identify clearly the 
contributions to E an i S resulting from level-shifts and lifting of degeneracies near Ep. In 
particular, the latter give important contibutions to E anis . We find for these that generally 
Eanis oc A^, is valid. The degeneracies do not only appear at high symmetry points of 
the Brillouin zone (BZ) |2j but also along lines in k-space. It is very important that we 
obtain convergent results for E an i S without applying state-tracking or excluding surface-pair 
coupling as was necessary in previous analyses by Wang et al. pj. Furthermore, we find the 
characteristic scale for the temperature dependence of the magnetic anisotropy to be A so , 
rather than the bandwidth, again indicating the significance of the lifting of degeneracies at 
Ep by A so . 

We start our theory by defining the magnetic anisotropy energy as 

E a nis{n) := E tot (6 = 0;n) - E tot (0 = 7r/2,0 o ;n), 

where E to t(6 , 4>] n) is the ground-state energy per atom of the monolayer with a total of n 
3d- and 4s-electrons per atom and A" atoms and is given by 

E tot (6, <P;n) = ^J2 ^mk(0, 4>)fo (Mn) - #mk(0, </>)) • 
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The angles 9 and <p denote the direction of the magnetization ( |§. fo(AE) is the Fermi- 
function at zero temperature and Ep{n) is the Fermi-energy, which, for a given bandfilling 
n, is determined self-consistently by 

m,k 

E m k(6, 4>) is the ra-th eigenvalue with crystal momentum k and magnetization along (9, 4>) 
of the Hamiltonian 

H = H s + H d + iJ srf + if so . 

H s refers to the s-electrons which are described by a set of properly symmetrized plane 
waves. Hd refers to the (^-electrons treated by spin-polarized tight-binding approximation 
using the Fletcher- Wohlfahrt parametrization |7|,§] adapted to the monolayer. H s d denotes 
s-d hybridization (in the conventional form of the combined interpolation scheme |J) and 
H s0 = A so l-s is the SOC between the (i-electrons on the same site in the usual approximation 
with the SOC parameter X so . H so is a matrix function of the magnetization direction (. It 
can be expressed in terms of the orbital momentum operators 1%, l v and with respect to 
the rotated frame (£, fj, Q as P,pH 
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For the numerical analysis, we employ for the ci-electrons a basis of 10 two-dimensional 
Bloch-functions ipi, i = 1, 10, for each possible crystal momentum k, constructed from 
atomic 3d-wave functions together with the spin eigenstates ||) and ||) of the Pauli matrix 
cr^, where ( is the spin quantization axis. Only nearest neighbors (4 in a simple- quadratic 
monolayer) are considered. To obtain accurate parameters, we perform a fit to the full- 



potential LMTO calculation for a free-standing Fe-monolayer by Pustogowa et al. |TT[ and 
to the LAPW calculation for a Ni- monolayer by Jepsen et al. [[HJ. The s- and d- bandwidths 
and hybridization parameters are then scaled to take into account the Cu surface lattice 
constant of 2.56 A ||13|| . The complete BZ summation over k is performed as a weighted 
summation over the irreducible part of the BZ (1/8 of the quadratic BZ in the nonmagnetic 
case or for perpendicular magnetization and 1/4 for in-plane magnetization due to reduced 



symmetry). To achieve convergence, about 1000 or 2000 points in the irreducible part of the 
BZ are sufficient. Note, we do not have to exclude any parts of the BZ to obtain convergence, 
unlike Wang et al. 0. 

In Fig. [l] we present results for E an i s as a function of the bandfilling n, in order to demon- 
strate the correspondence between electronic structure and magnetic anistotropy energy and 
to show that our method will yield convergent results for the whole transition metal series 
and for large (Fe) and small (Ni) exchange coupling. We use parameters for freestand- 
ing Fe and Ni monolayers and a lattice constant of 2.56 A to simulate epitaxial growth 
on Cu(001). In particular, our results yield that Fe has a perpendicular easy axis with 
Eanis — —0.32 meV and Ni an in-plane easy axis (along an axis connecting nearest neigh- 
bors) with E an i S = 0.1 meV ||14| . Experiment yields perpendicular anisotropy for ultrathin 
Fe-films 0. However, comparison to experiment for a monolayer is difficult due to growth 
problems. In the case of Ni, even the calculated in-plane direction of the easy axis agrees 
with experiment M. Note, corresponding ab initio results for a free-standing Fe-monolayer 
yielded —0.4 meV IPUT5I, but previous tight-binding calculations gave —5.5 meV fll6| . Al- 
though our numerical values of E anis depend on the choice of parameters, the sign and the 
order of magnitude of E an i S are remarkably stable upon parameter variations: In agree- 
ment with Wang et al. [[jj we find a perpendicular easy axis also for Fe monolayers taking 
(001) surface lattice constants of Pd, Ag, V, and W (2.77 A, 2.89 A, 3.03 A, and 3.16 A, 
respectively) . 

It is of considerable interest for the physical interpretation of the origin of magnetocrys- 
talline anisotropy energy to study E anis as a function of A so and the hopping parameters t. 
To a good approximation, we find quite generally E anis oc Af , which sheds light on previous 
analyses. This was calculated for Fe and Ni bandstructure parameters for fixed bandfilling n 
and in the range of A so from zero to the actual value (Fe: A so = 50 meV, Ni: A so = 70 meV). 
However, for very small E an is the dependence might be larger that \ 2 so . 

Regarding the dependence of E an is on A so , the following remarks are of interest. Treating 
SOC as a perturbation, the lowest nonvanishing order of non-degenerate perturbation theory 
is \ 2 so due to time reversal symmetry |T0| . However, in the case of degenerate bands, a lifting 
of the degeneracy by SOC will complicate the situation. Degenerate bands can undergo a 



4 



lifting linear in \ so when SOC is introduced, and the resulting contibution to E anis depends 
on the area in k-space influenced by the degeneracy. Whether this area is of the order of 
\ 2 so , which would yield E anis oc A^ |J, or this area is of lower order and thus would yield an 
important contribution to E an i S [@,[n]L has been a controversial question. This makes our 
result very interesting (see below). 

Concerning the dependence of E anis on the (i-electron hopping parameters t, note that 
the overall shape of the curves E anis (n) will not change if t is varied. | .Earns I increases 
for decreasing t (decreasing bandwidth). This leads to the general trend: \E anis \ increases 
with increasing lattice constant of the monolayer. This corresponds to an increasing sub- 
strate lattice constant, if one neglects hybridization effects between the substrate and the 
monolayer. 

For the further physical interpretation of our results for E an i S as a function of n and 
X so , namely how E anis results from the bandstructure and especially how it is affected by 
the lifting of degeneracies in the bandstructure close to Ep, we neglect s-d hybridization 
for simplicity. This is allowed since only d-state degeneracies matter; s-d hybridization 
changes only the location of the degeneracies. The inset in Fig. |2] shows schematically that 
a special type of degeneracy ( "line" degeneracy) leads to an important contribution to E anis 
if it occurs near Ep for one direction of magnetization and is lifted by SOC for another. 
Estimating this contribution (see Fig. [|) yields AE an i S = |A so -F. The fraction F of involved 
states in k-space is given by F = ^r 1 • 1 = 2X so (f fr-) , since the intersecting bands are 

a \ 1 / 

non-dispersive along one direction in k-space (perpendicular to ki). Here, a is the lattice 
constant, and is the dispersion of the intersecting bands near their intersection. Hence, 
the maximum energy gained by the lifting of the degeneracy amounts to 

It is important to note that such contributions are not caused by degeneracies at isolated 
points [I8|, but by degeneracies along lines in k-space. Our bandstructure calculations yield 
such degeneracies near Ep for n = 8.7 using Fe parameters for the bandstructure and for 
n = 5.9 using Ni parameters. 

Eq. (0) explains immediately the height of the peak at n = 8.7 shown in Fig. [I] (curve 
a). Taking from our bandstructure calculations (Fe parameters) ~ 1 eV/^ we obtain 



AE an i S m 2.5 meV. This is in excellent agreement with the exact result of our calculation 
E an i S (n = 8.7) = 2.7 meV. In addition, a detailed k-space resolved calculation of E anis {n = 
8.7) shows that indeed only states near the corresponding degeneracy at Ep yield large 
contributions to E an i S . Furthermore, Eq. (FJ) shows that AE an i S oc \ since approximately 
^ oc t. Note, this is in agreement with our previous observation regarding the dependence 
of E anis on the hopping parameters t. 

Since the degeneracies are so important for E anis , we outline in the following at some 
length the relationship between E an i S and the degeneracies and thus the crystal symmetry. 
For that purpose, we adapt our basis set of Bloch wave functions tpi, i = 1, 10, in order to 
obtain a maximum number of zeros in the Hamiltonian matrix for the (i-electrons. Expressed 
in terms of Bloch functions constructed from the cartesian atomic orbitals, the 5x5 tight- 
binding matrix has its simplest block diagonal form with only two off-diagonal elements 
(ODEs) (-Hrf)4,5 = (-Hd)5,4 (and equivalently (H d ) 9 i = (H d ) 10 $) if the x-axis is directed along 
an axis connecting nearest neighbors. To find out which additional ODEs are introduced by 
SOC for a given direction of the magnetizaton M (in the following, M || z and M || x are 
considered), we analyze the form of H so in Eq. (JJ). We assume that different Bloch-states 
are orthogonal. States with parallel spins are coupled, if they contain equal orbital momenta 
with respect to the spin quantization axis (, whereas states with opposite spins must show a 
difference of one in the orbital momenta to yield nonvanishing ODEs. The cartesian orbitals 
xy, yz, zx, x 2 — y 2 and 3z 2 — r 2 (leading to ipi, ...,"05 together with |f) and to ip e , ...,ipio 
with ||), respectively) are composed of eigenstates of l z with the eigenvalues (-2,2), (-1,1), 
(-1,1), (-2,2) and 0, respectively. In terms of eigenstates of l x one finds the eigenvalues (-1,1), 
(-2,2), (-1,1), (-2,0,2) and (-2,0,2), respectively. Hence, coupling exists, for M || z, within 
the groups of states ^ with i = 1,4,5,7,8 and with i = 2,3,6,9,10, and, in the case of 
M || x, within the groups of states ipi with i = 2,4, 5, 6, 8 and i = 1,3, 7, 9, 10, respectively. 
In both cases, the Hamiltonian can be split into two 5x5 blocks, and subbands belonging 
to different blocks will intersect. Between states of the same block, the degeneracies will 
ordinarily be removed. Especially the subbands ip\ and ip2 (and, correspondingly, ip G and 
ipi) change their roles, if the magnetization is changed from z to x and vice versa, because 
the orbitals xy and yz have different orbital momenta with respect to the x- and z-ax.es. So, 
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they will be involved in the lifting of degeneracies by altering magnetization and possibly, as 
shown above, yield important contributions to E anis . As an example, the prominent peak in 
the E anis (n) curve of Fe at n = 8.7 results from a degeneracy of the subbands corresponding 
to the states ^7 and (^9, tpio) for M || z, occurring along a line parallel to the k x -&xis in 
k-space, that is lifted for M || x, because in the second case the subbands belong to the 
same block of the Hamiltonian, whereas in the first they do not. 

Concerning the different behavior of E anis as a function of n choosing Fe or Ni parameters 
(see Fig. |I]), one should note that due to the large exchange splitting in Fe only the SOC 
between d-states with parallel spins contributes significantly to E anis . Thus, E anis for a less 
than half-filled <i-band (n < 6.2 in Fig. |I|, curve a) mainly results from the majority-spin 
subband and E an i S for a more than half-filled d-bernd (n > 6.2) from the minority-spin 
states. This explains the similar shape of these two parts of the curve E an i S (n). For Ni, the 
contribution to E anis from SOC between opposite spin states is of equal magnitude, leading 
to a rather different behavior of E anis (n). 

Apparently, the occurrence and lifting of degenerate subbands for different directions of 
magnetization is important for calculating the magnetic anisotropy energy. This should also 
be clearly seen from the temperature dependence of E anis due to the magnetization M(T), 
the Fermi function Jt{AE), the hopping integrals t(T), which depend on temperature due to 
lattice expansion, and the entropy S(T). In particular, effects resulting from the latter three 
contributions are analyzed. We find a characteristic energy scale of the order of magnitude 
of X so for the reduction of the free anisotropy energy F anis = E anis — TS anis with increasing 
temperature. This is demonstrated in Fig. ^| in a <i-band calculation with Fe parameters 
(d-bandfilling = 6.0), where F anis decreases over a temperature range of 500-1000 K 
(50-100 meV), which corresponds to the energy X so , but not to the 3<i-bandwidth of about 
3 eV. This becomes immediately plausible if one notices that the SOC-induced lifting of 
degeneracies occurs near the Fermi-level. Thus one expects a measurable effect on F an i s if 
/cbT becomes larger than 2A so . Furthermore we must conclude from our results that shifting 
of subbands far below the Fermi-level is not so important, since then F anis could not be 
essentially lowered on such a small temperature scale. Hence, this analysis of F an i S (T) shows 
also the significant role of changes of the degeneracies giving an important contribution to 
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the anisotropy energy. It is remarkable that the three temperature effects mentioned above 
are of equal magnitude as the temperature effects due to M(T). 

In conclusion, a fully convergent calculation of the magnetocrystalline anisotropy energy 
Eanis of Fe and Ni monolayers on Cu (001) is performed. We find a perpendicular easy axis 
for Fe and an in-plane easy axis for Ni. Large contributions to E anis result from the SOC- 
induced lifting of degeneracies at the Fermi-level, which is also supported by our calculation 
of temperature effects. In general, E anis scales with the square of the SOC constant A so . 
The main features of our analysis can be applied also to thicker films, which is of particular 
interest for the investigation of reorientation transitions. 
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FIGURES 

FIG. 1. Dependence of E an i s on the 3d- and 4s— bandfilling n for a monolayer with parameters 
referring to Fe (curve a) and Ni (curve b). Negative values of E an i s yield perpendicular anisotropy. 
The vertical lines denote n for Fe and Ni, respectively. 

FIG. 2. Temperature dependence of F an i s (T) for a Fe-parametrized ci-band calculation with 
d-bandfilling rid = 6. The inset shows the occurrence (full lines) and lifting (dashed lines) of 
a "line" degeneracy for two different directions of magnetization Mi and M2, respectively, ki 
corresponds to one particular direction in k-space. Perpendicular to ki the intersecting bands are 
non-dispersive throughout the BZ. Note, the energy gained by the lifting of this degeneracy is given 
by AE an i s = ^X so -F, if Ep falls in between the two subbands (dotted line). Here, F is the fraction 
of the involved states in k-space. Apparently, if Ep lies below or above the two subbands, AE an i s 
is zero. 
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